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Abstract 

An engineering inviscid-boundary layer method 
has been modified for application to slender three- 
dimensional (3-D) forebodies which are characteristic of 
transatmospheric vehicles. An improved shock descrip- 
tion in the nose region has been added to the inviscid 
technique which allows the calculation of a wider range 
of body geometries. The modified engineering method is 
applied to the perfect gas solution over a slender 3-D con- 
figuration at angle of attack. The method predicts sur- 
face pressures and laminar heating rates on the windward 
side of the vehicle that compare favorably with numeri- 
cal solutions of the thin-layer Navier-Stokes equations. 
These improvements extend the 3-D capabilities of the 
engineering method and significantly increase its design 
applications. 

Nomenclature 

A, B, C, D, E circumferential conic section parameters 
b, c, d longitudinal conic section parameters 

e,f,g shock curvature parameters 

e^,ey unit vectors of shock-oriented coordinate 

system 

L reference length 

M Mach number 

n coordinate normal to shock 

p pressure 

q heat transfer rate 

s distance 

V velocity magnitude 

V velocity vector 

x,r,<f) cylindrical coordinate system 

x,y, z Cartesian coordinate system 

a angle of attack 
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k shock curvature 

v shock angle in a meridional plane 

£,/? shock coordinates 

p density 

stream function 

Subscripts 

b body 

i marching step 

k meridional plane 

sub subsonic region 

w wall 

oo freestream conditions 


Introduction 

The aerothermal design of hypersonic vehicles con- 
sists of evaluating vehicle geometries for a range of flow 
conditions. Numerical solutions of the Navier-Stokes 
(NS), parabolized Navier-Stokes (PNS), or viscous shock- 
layer (VSL) equations are usually cost-prohibitive except 
for the later stages of the design. For preliminary design 
and parametric studies, engineering inviscid-boundary 
layer methods are more appropriate. These methods pre- 
dict reasonably accurate surface heating rates at a frac- 
tion of the computer cost of the NS algorithms. 

An improved engineering inviscid-boundary layer 
method 1, 2 has been developed by combining a simpli- 
fied three-dimensional (3-D) inviscid technique 3, 4 with 
the axisymmetric analog concept 5 and a set of approx- 
imate convective heating equations. 6 This algorithm is 
referred to as THINBL (Three-dimensional Hypersonic 
INviscid-Boundary Layer method). The inviscid tech- 
nique, which is similar to the thin shock-layer method 
of Maslen, 7, 8 improves the surface pressure distribution 
and extends the 3-D capabilities of existing engineering 
methods. The method is applicable to perfect gas and 
equilibrium air flows at both laminar and turbulent con- 
ditions. Solutions have been obtained for axisymmetric 
shapes (paraboloids, sphere-cones) and simple 3-D shapes 
(blunted elliptic cones) at angle of attack. 1 Comparisons 
with experimental data and NS solutions have shown 
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Figure 1. Slender 3-D forebody of a transatmospheric 
vehicle. 

good agreement. This engineering method provides a 
unique capability for approximate flowHeld solutions over 
3-D blunted noses as well as 3-D afterbodies. 1 

The next step in the development of the THINBL 
code is to apply it to a realistic vehicle geometry. One 
such geometry is the slender 3-D forebody characteristic 
of a transatmospheric vehicle (see Fig. 1). However, the 
first attempt at computing the flowHeld over this body 
highlighted several limitations of the algorithm. The 3-D 
inviscid technique used by THINBL is an inverse method: 
an assumed shock shape is iterated until t.hq Actual body 
shape is calculated. In the subsonic region near the nose, 
the possible shock surfaces assumed by THINBL cannot 
produce the slender 3-D forebody shown in Fig. 1. For 
application to this configuration, the shock description 
requires modification. In addition, the sonic line loca- 
tion (which determines the subsonic boundary and the 
location win re I In solution can be marched downstream) 
varies greatly around the circumference of the body. This 
variation is currently neglected in the 3-D inviscid tech- 
nique. Also, due to the large curvature near the sides of 
the vehicle, some means of redistributing streamlines is 
necessary. 

The purpose of this paper is to address these lim- 
itations and extend the capabilities of the engineering 
method THINBL to include more complex geometries 
such as the slender 3-D forebody of a transatmospheric 
vehicle. This would significantly increase the design ap- 
plications of engineering methods. Surface pressures and 
heating rates are computed for this geometry at angle 
of attack. Comparisons are made with NS solutions to 
demonstrate the accuracy of the present technique and 
to document t.h# modifications. 

Analysis 

This .section describes the engineering inviscid- 
boundary layer method THINBL. The modifications 


Figure 2. Shock-oriented curvilinear coordinate system: 
side view. 

needed for the solution over slender hypersonic config- 
urations are also discussed. 

THINBL 

The engineering code THINBL (Three-dimensional 
Hypersonic INviscid-Boundary Layer method) is com- 
prised of an approximate 3-D inviscid method, a proce- 
dure for computing inviscid surface streamlines, and a 
set of approximate convective heating equations. The ax- 
isymmetric analog concept is employed for calculating the 
heating rates along surface streamlines. A complete de- 
scription of THINBL is given in Refs. 1 and 2. A brief 
outline is presented here. 

Inviscid Method 

The 3-D analysis is simplified by writing the inviscid 
equations in a shock-oriented curvilinear coordinate sys- 
tem as shown in Figs. 2 and 3 where ^ and f3 represent 
coordinates of a point on the shock surface and n (not 
shown) is the inward distance normal to the shock. The 
unit vectors, and e^, are tangent to the shock surface 
and are chosen such that is in the direction of the tan- 
gential velocity just inside the shock surface. The unit 
vector e /3 is then defined to be perpendicular to and 
e n . A simplifying assumption is made that the velocity 
component in the direction is equal to zero not only 
at the shock but throughout the shock layer. This ap- 
proximation is similar to the axisymmetric analog used 
in the boundary-layer analysis and allows the solution 
along a given £-line to be computed independent of other 
lines. The momentum equations in the curvilinear coor- 
dinate system are then simplified by following the thin 
shock layer analysis of Maslen. 1 ’ 8 This yields explicit ex- 
pressions for the pressure and normal velocity component 
along the n coordinate as a function of a stream function 
'F. 
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Figure 3. Shock-oriented curvilinear coordinate system: 
front view. 

As stated previously, this technique is an inverse 
method in which the shock is assumed known. Flowfield 
properties are computed along n from the known shock 
(\? = \? s ) to \? = 0 which defines the body. The shock 
shape is iterated until the computed body shape matches 
the actual body geometry. In the subsonic-transonic re-; 
gion, the shock is globally iterated until the correct body 
shape is obtained. Beyond this region, the shock surface 
is locally iterated as the solution is marched downstream 
along a given number of ,f-li nes - 

Axisymmetric Analog 

The 3-D boundary-layer analysis is simplified by us- 
ing the axisymmetric analog 5 as is done in most engineer- 
ing aerothermal methods. The 3-D boundary-layer equa- 
tions are first, written in a inviscid surface streamline co- 
ordinate system. The crossflow velocity component tan- 
gent to the surface but normal to the streamline is then 
assumed to be equal to zero throughout the boundary 
layer. This simplification reduces the 3-D boundary-layer 
equations to the axisymmetric form provided the distance 
along the streamline is substituted for the surface dis- 
tance and the scale factor describing the divergence of 
the streamlines is interpreted as the axisymmetric body 
radius. Axisymmetric boundary-layer methods can then 
be coupled to the 3-D inviscid technique. 

Before applying the axisymmetric analog, inviscid 
>n r face streamlines are computed from the approximate 
inviscid solution. These body streamlines may be cal- 
culated from the surface pressure distribution 9 or from 
the velocity components. 10, 13 The approximate inviscid 
method used in THINBL predicts accurate surface pres- 
sures, but the direction of the velocity on the surface is 
not accurate due to the assumption of zero velocity in 
the ep direction. Therefore,; streamlines are calculated 
from the surface pressures in THINBL. The equations for 
computing th§ streamlines are given in Refs. 1 and 2. 


Although the surface streamlines can be determined af- 
ter the inviscid solution has already been calculated, it is 
more convenient to compute the inviscid solution and the 
streamlines simultaneously. In THINBL, streamline in- 
formation (such as its orientation and the scale factor de- 
scribing the divergence of streamlines ) is integrated along 
a £-line along with the inviscid flowfield variables. 

Boundary-Layer Method 

The axisymmetric analog allows any axisymmetric 
boundary-layer method to be applied along an inviscid 
surface streamline. In THINBL, a set of approximate 
convective-heating equations developed by Zoby b is used 
for the boundary-layer solution. Laminar and turbulent 
heating rates may be calculated from these relations for 
both perfect gas and equilibrium-air flows. Approximate 
expressions for the boundary-layer thickness at both lam- 
inar and turbulent conditions are also given in Ref. 6. 
Boundary-layer edge conditions are found by interpolat- 
ing in the approximate inviscid solution a distance from 
the wall equal to the boundary-layer thickness. This ap- 
proach has been demonstrated to approximately account 
for the effects of entropy-layer swallowing. 

Modifications 

To compute the flow field over the slender configu- 
ration in Fig. 1, several areas of THINBL require modi- 
fication. Changes in the shock shape near the nose; the 
handling of the subsonic region; and the distribution of 
the teli nes around the body are discussed. 

Shock Shape 

In the subsonic region surrounding the nose of the 
vehicle, the shock surface must be globally iterated. 
Presently, the shock surface is described by three longi- 
tudinal conic sections with the circumferential variation 
given by an ellipse. In a meridional plane, the longitudi- 
nal fit is given by the conic section 

( r 2 + bx 2 — 2cx + ‘2dxr)k = 0 (1) 

where k represents the local meridional plane and the 
parameters (b,c,d.)k govern the shape of the conic. 
The cylindrical coordinates (x,<f>,r) are defined in a 
wind-oriented system with the x axis aligned with the 
freest.ream velocity. Since one ellipse is used for the cir- 
cumferential variation of the shock, only three meridional 
planes (<ft = 0,90, 180 deg) are needed to completely de- 
scribe the shock. The ellipse is given by 

r 2 ( A cos 2 <p + sin 2 c/>) + rC cost i> = E (2) 

white A, G, E depend on the local values of rj, from the 
three meridional planes and are functions of x. This form 
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of the shock produces a body that can similarly be de- 
scribed by conic sections with an elliptical cross-section. 
However, the cross-section of the slender 3-D forebody 
shown in Fig. 1 cannot be described by a single ellipse 
nor can the body be accurately fitted with a conic sec- 
tion in the planform view. Changes are required in both 
the axial and the cross-sectional fits of the shock surface 
for application to this vehicle. 

A more general axial fit for the shock is the specifica- 
tion of the longitudinal shock curvature with a quadratic 
in x as 

K k = (e + fx + gx 2 ) k (3) 

where k is the curvature. The radius is computed by an 
integration of 


d sin v k 
dx 


= ~K k 


dr k 

—— = tan v k 
dx 


( 4 ) 


where taniz represents the shock slope. This fit was sug- 
gested by Maslen 11 and produces conic as well as non- 
conic body shapes. The cross-sectional fit is improved by 
matching conic arcs. Each conic arc is described by 


r 2 (Acos 2 <f> + B sine/) cos (/>+ sin 2 (/>) + r(C cos^ + D sin^) = E 

( 5 ) 

At the match points between the arcs, the radius and its 
first and second derivatives with respect to <f> are forced 
to be continuous. Cheatwood and DeJarnette 12 follow a 
similar approach in curve-fitting body geometries except 
that no constraint is placed on the second derivative of the 
radius. The addition of conic arcs increases the number 
of parameters that describe the shock and allows for the 
calculation of body geometries with non-elliptical cross- 
sections. 


Subsonic Region 

The global iteration of the shock surface in the sub- 
sonic region involves matching the body shape computed 
from the approximate flowHeld equations to the actual 
body geometry. To insure a good starting solution for 
the downstream marching procedure, the matched body 
points are located at the end of the subsonic region where 
the flow on all £-lines is supersonic. For simplicity, the end 
of the subsonic region in the THINBL code is assumed to 
be located at x = x su b where x su b is constant and en- 
compasses the entire subsonic region. The variation of 
x S ui with <f> is neglected for the previous applications to 
a blunted elliptic cone at angle of attack. 

However, the end of the subsonic region is not as con- 
veniently defined for the case of a slender 3-D forebody. 
The flow near the “edges” of the configuration remains 
subsonic at much longer distances than the flow near the 
lower and upper symmetry planes. It is impractical to 
globally iterate the shock surface over this entire region 
with only a few parameters. For this reason, x su b is now 
allowed to vary for each of the meridional planes used in 
the shock fit. 


^-Line Distribution 

One of the difficulties in tracing £-lines is that they 
tend to wrap around the 3-D shock surface. A uniform 
distribution of £-lines is not guaranteed. For the previous 
applications to axisymmetric shapes and blunted ellip- 
tic cones at angle of attack, this problem is handled by 
grouping the £-lines in the stagnation region such that an 
even distribution is obtained at some point downstream. 
However, for a slender 3-D forebody, the flow is highly 
three-dimensional near the edges of the configuration and 
£-lines diverge rapidly from this region. 

To insure a uniform distribution of £-lines as the so- 
lution is marched downstream, £-lines are redistributed 
at every axial station. Inviscid flowHeld and streamline 
variables are integrated along a number of lines from Xi 
to A new set of £-lines is chosen at X {. |_i that uni- 

formly encompasses the shock. FlowHeld and streamline 
variables associated with the new set of £-lines are found 
by interpolation between the previous lines. The process 
is repeated at each subsequent marching step. This ap- 
proach has been suggested and used by Hamilton 13 to 
compute the streamlines and surface heating over a com- 
plex vehicle. 


Results and Discussion 


Surface pressure and heating rate distributions are 
presented at laminar, perfect gas conditions over an ax- 
isymmetric power-law body and a slender 3-D forebody 
at angle of attack to assess the improvements and accu- 
racy of the THINBL code. All variables are nondimen- 
sionalized as follows: pressure by PcoV^, heating rate by 
PcoV^, and distances by a reference length L. 

Axisymmetric Power-Law Body 

To examine the accuracy of the longitudinal shock 
fit for non-conic body shapes, the flowHeld is computed 
over an axisymmetric power-law geometry given by 

n = x\ /A ( 6 ) 

This geometry models the planform view of the slender 
3-D configuration in Fig. 1, and the configuration can- 
not be described adequately by a conic section. The 
freestream conditions are = 15, poo = 0.00949 

kg/m 3 , and Tx, = 234 K. The wall temperature is 
T w = 1256 K and the reference length is L = 0.10 m. Fig- 
ure 4 gives the body shape calculated by the conic section 
shock, and Fig. 5 shows the body geometry which is pre- 
dicted by assuming that the shock curvature is given by 
a quadratic in x. Excellent agreement between the com- 
puted and actual body geometries is given by the modified 
shock shape. The shock curvature predicted by both lon- 
gitudinal fits is presented in Fig. 6. The quadratic fit al- 
lows for a much smoother transition between the subsonic 
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Figure 7. Surface pressure comparison for axisymmetric 
Figure 4. Comparison of body shapes produced by conic power-law body 
section shock. 
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Figure 5. Comparison of body shapes produced by 
quadratic for shock curvature. 
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Figure 6. Shock curvature distribution. 


region (x < 0.25) and the supersonic region downstream. 

Nondimensional surface pressures and heating rates 
are shown in Figs. 7 and 8 for the axisymmetric power- 
law geometry. The distance s b is the distance from the 
origin and is given by 

■s h = [x 2 b +rj} 1/2 (7) 

Near the nose, s b approximates the distance along the 
surface. Results of THINBL with the quadratic and 
conic shock shapes are compared with a solution given 
by the thin layer NS algorithm LAURA. 14 Because of 
the error between the calculated and actual body shapes, 
THINBL with the conic shock shape overpredicts the 
pressures generated by LAURA by as much as 35 per- 
cent. The modified version that uses the quadratic for 
the shock curvature gives more accurate (within 10 per- 
cent) and smoother pressures in the nose region. Note 
that both versions predict the same pressures downstream 
(s b > 0.75). In Fig. 8, the surface heating rates predicted 
by the modified THINBL follow the same trends as the 
heating rates from LAURA. The heating increases away 
from the stagnation point due to the shock layer becom- 
ing thinner. The original version of THINBL does not 
account for this effect. 


Slender 3-D Forebody 

To examine the accuracy of the improved 3-D shock 
description, the flowHeld is computed over a slender 3-D 
forebody at 10 deg angle of attack. This geometry is 
characteristic of a transatmospheric vehicle. The side 
view is described by a circular nose followed by a line in- 
clined at 5 deg; the planform is described by a power law; 
and the cross-section is given by a super-ellipse where 
the cos 2 <f) and sin 2 (/> terms in Eq. 2 have been replaced 
by cos 4 (/> and sin 4 (/>. All three views of i In; vehicle are 
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Figure 8. Surface heating comparison for axisymmetric 
power-law body. 



Figure 9. Slender 3-D forebody: side view. 

shown in Figs. 9-11. Hole that the three views are not 
drawn to the same scale. The freestream conditions are 
Mco = 10, pco = 0.00170 kg/m 3 , and XA = 48.3 I\. The 
wall temperature is T w = 288 I\ and the reference length 
is L = 1.905aT0 -4 m. Surface pressures and heating rates 
generated by THINBL are compared with thin layer NS 
results from the LAURA algorithm. Ninety-one ^-lines 
are used in THINBL. They are distributed only on the 
lower half of the vehicle because of the code’s inability to 
handle the separation present on the leeside. The ^-lines 
are redistributed at each axial location to control their 
divergence from the sides of the vehicle. A grid of 40 
cells in the axial direction, 200 cells around the circum- 
ference of the body, and 64 cells in the normal direction 
is used to obtain the LAURA solution. Both the ^-lines 
in THINBL and the circumferential cells in LAURA are 
clustered in the high curvature regions near t/j = 0 (See 
Fig. 11). THINBL requires 15 minutes to obtain a so- 
lution on a Sun Sparcstation while the LAURA solution 
required 25 hours on a Cray Y-MP supercomputer. 



Figure 10. Slender 3-D forebody: top view. 



Figure 11. Slender 3-D forebody: rear view. 

Axial surface pressures and heating rates are shown 
in Figs. 12 and 13 for the windward plane of symmetry. 
The Mirface pressures compare well with the pressures 
from the LAURA solution on the nose, but underpredict 
the pressures aft, of the nose by 25 percent. Similar dis- 
crepancies have been seen downstream of the sphere-cone 
juncture of blunted slender cones. 4 One of the thin shock- 
layer assumptions in THINBL ’s inviscid technique is that 
the local shock and body slopes are approximately equal. 
In the region near the nose of the slender 3-D forebody 
and near the sphere-cone juncture of a blunted slender 
cone, this assumption is inaccurate. As a result, the sur- 
face pressures are in error in this region. However, for 
a blunted slender cone, the errors diminish after a few 
nose radii as shown for a 5 deg blunted cone in Fig. 14. 
Over the slender 3-D forebody, the differences between 
the THINBL and LAURA pressures persist for the 15 
nose radii presented in Fig. 12. Because the lower surface 
of this 3-D configuration is flat, the centerline flowheld is 
almost two-dimensional (2-D). The shock layer is thicker 
in a 2-D flow than in an axisymmetric flow. As a result, 
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Figure 12. Axial pressure comparison along windward 
centerline of slender 3-D forebody. 



Figure 13. Axial heat, transfer comparison along wind- 
ward centerline of slender 3-D forebody. 

the accuracy of the thin shock-layer assumptions is worse 
for a blunted slender wedge than for a cone. However, 
the approximations improve downstream as depicted in 
Fig. 15 for the inviscid flow over a 5 deg blunted wedge. 

It is felt that the comparison of centerline pressures over 
the slender 3-D vehicle will improve downstream as well. P w 

Axial surface heating rates given by THINBL com- 
pare favorably (within 15 percent) with the LAURA so- 
lution in Fig. 13. Like the pressure comparison, the heat- 
ing rates from THINBL underpredict the LAURA results 
downstream of the nose. These comparisons should again 
improve farther downstream. Preliminary comparisons 
with LAURA pressures and heating rates in the wind- 
ward symmetry plane farther downstream confirm this 
hypothesis. 

Circumferential pressures and heating rates are pre- 
sented in Figs. 16 and 17 for an axial location of x = 10. 

The surface variables are shown as a function of rj where 
zi, = 0 represents the upper or lower symmetry plane 



Figure 14. Axial pressure comparison for 5 deg blunted 
cone. 



Figure 15. Axial pressure comparison for 5 deg blunted 
wedge. 



Figure 16. Circumferential pressure comparison for 
slender 3-D forebody. 
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Figure 17. Circumferential heat transfer comparison for 
slender 3-D forebody. 



Figure 18. Circumferential pressure comparison for 
slender 3-D forebody. 

and zi, progresses toward the outermost point on the ve- 
hicle surface. The pressures and heating rates are essen- 
tially constant on the lower surface of the vehicle and 
increase dramatically as the edge is approached. Good 
agreement between the THINBL and LAURA solutions 
is seen on the flat lower surface and near the outer edge of 
the vehicle. A closer examination of the properties near 
the vehicle’s edge is given in Figs. 18 and 19. The pres- 
sures and heating rates are presented with respect to (f>b 
with (f>b = 90 representing the edge of the vehicle. Good 
agreement (within 15 percent) between the THINBL and 
LAURA pressures is seen in this region. An interesting 
feature is noted in the heating rate comparisons in Fig. 19. 
Both THINBL and LAURA predict two peaks in the heat- 
ing in this area. The peaks probably result from the fact 
that the super-ellipse cross-section of the vehicle is flat, 
at <f>b = 90. Thes; maximum curvature of the body cross- 
section is found at (f>b < 90 and (f>b > 90 (see Fig. 11). The 
surface streamlines diverge rapidly from these regions of 



Figure 19. Circumferential heat transfer comparison for 
slender 3-D forebody. 

high curvature. This streamline divergence may locally 
increase the heating levels. The magnitude of the heating 
rates from THINBL compare favorably with the LAURA 
solution for the first, (and largest.) peak, but. underpredict, 
the heating just, on the leeside (<pb > 90). This ability to 
compute flowfield information and surface heating rates 
a.wa.y from the plane of symmetry for a. 3-D vehicle in- 
creases THINBL ’s applicability as a. preliminary design 
tool. 
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Concluding Remarks 

An engineering inviscid-bounda.ry layer method has 
been modified for application to hypersonic, slender 3-D 
vehicle forebodies. A new description of the shock shape 
in the nose region has been included. The engineering 
method with the new shock is applied to the solution over 
an a.xisymmet.ric power-law shape and a. slender 3-D con- 
figuration at. angle of attack. Comparisons with thin layer 
NS solutions demonstrate the accuracy of the new shock 
shape and its improvement, over the original version. Sur- 
face pressures and heating rates compare favorably with 
the NS solution on the centerline and on the lower surface 
of the vehicle. The underpredict. ion of the pressures and 
heating rates on the centerline aft. of the nose should im- 
prove farther downstream. These modifications represent, 
a. significant, improvement, over current, engineering meth- 
ods and enhance current, approximate 3-D capabilities for 
realistic slender geometries. 
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